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We investigate vortex structure in very anisotropic layered superconductors at fields tilted with 
respect to the c-axis. We show that even a small in-plane field does not tilt the vortex lattice but, 
instead, penetrates inside the superconductor in the form of Josephson vortices (JVs). At high c-axis 
magnetic field the phase field of the JV is built up from the phase perturbations created by displace- 
ments of pancake vortices. The crossing-lattices ground state leads to linear dependencies of the 
melting field and melting temperature on the in-plane field, in agreement with recent experimental 
observations. At small fields stacks of JVs accumulate additional pancake stacks creating vortex 
rows with enhanced density. This mechanism explains the mixed chains-lattice state observed by 
Bitter decorations. 



Significant progress has been achieved recently in un- 
derstanding of vortex states and melting transition in lay- 
ered superconductors with very weak Josephson coupling 
such as Bi 2 Sr2CaCu20 a; (BSCCO). A magnetic field ap- 
plied perpendicular to the layers penetrates such super- 
conductors in the form of pancake vortices (PVs) Q| . P Vs 
in different layers are coupled weakly via the Josephson 
and magnetic interactions and form aligned stacks at low 
fields and temperatures (PV stacks). The magnetic na- 
ture of coupling leads to an unusual melting scenario. In 
contrast to superconductors with strong Josephson cou- 
pling, both crystalline order and alignment of PVs are 
destroyed at the melting point leading to almost inde- 
pendent two-dimensional vortex liquids in the layers Q . 

A rich variety of vortex structures were proposed for 
the magnetic field tilted with respect to the c-axis, such 
as the kinked lattice crossing lattices of Abrikosov 

and Josephson vortices (JVs) (combined lattice) j|, 
tilted vortex chains ||, etc. Obviously, the thermody- 
namics of the melting transition at tilted fields is deter- 
mined by the nature of the ground state. In moderately 
anisotropic superconductors an in-plane field simply tilts 
the vortex lattice inside the superconductor. This sce- 
nario does not work in the case of very weak Josephson 
coupling when the interlayer coupling is predominantly 
magnetic. The in-plane magnetic field interacts with the 
PVs only via the Josephson coupling. On the other hand, 
the PVs are mainly aligned by magnetic coupling, there- 
fore a homogeneous tilt of the lattice costs magnetic en- 
ergy. To save the magnetic coupling energy it is more 
favorable for the in-plane field to penetrate inside the su- 
perconductor in the form of JVs as in the Meissner state. 
The ground state of the system is then given by the di- 
lute lattice of the JVs coexisting with the dense lattice 
of the PVs (crossing lattices or combined lattice). This 
ground state was first proposed for the opposite case of 
a small c-axis field || when the dilute pancake lattice 



coexists with the dense Josephson lattice. The JV in 
weakly coupled superconductors has a very wide nonlin- 
ear core (the region where the interlayer phase difference 
is large). This means that even at moderate c-axis fields 
the JV core contains many PVs (see Fig. [l]a). We show 
that in this situation the phase field of the JV is built 
up from the phase perturbations created by pancake dis- 
placements. Such a JV has a smaller core size and smaller 
energy as compared to the ordinary JV. In the crossing- 
lattices state the free energy depends almost linearly on 
the in-plane field. As a consequence, the melting temper- 
ature and melting field of the pancake lattice also depend 
on the in-plane field in a linear way in agreement with 
recent experiments HQ. At high in-plane field the JV 
cores start to overlap. In the regime of strongly over- 
lapping JVs the in-plane field produces only weak zigzag 
deformations in the PV lattice Ejl. These deformations 
have a very little influence on the phase distribution and 
can be treated perturbatively, in contrast to isolated JV. 

At very small c-axis fields the JV core contains only 
one row of PVs. In this case displacements of the PVs 
have a very little influence on the JV structure. These 
displacements, however, produce a finite interaction en- 
ergy between the PV stack and the JV. Due to a very 
high anisotropy the JV lattice consists from densely 
packed stacks along c-axis (JV stacks) separated by rel- 
atively large distances. PV stacks crossing JV stacks 
have smaller energy as compared to other stacks and it 
would be favorable to add extra PV stacks to the rows 
located on the JV stacks. However this makes the PV 
lattice defective and increases its energy. On the other 
hand, when the distance between PVs a is much larger 
than the London penetration depth A, the energy of de- 
fect in the lattice is exponentially small. This means 
that at a certain field there is a phase transition from 
the weakly deformed triangular lattice to the defective 
lattice. In the latter state PV rows located on the JV 
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stacks will have a larger density of PVs than other rows 
(see Fig. |l]b). This picture provides a natural explana- 
tion for the mixed chain-lattice state observed by Bitter 
decorations in BSCCO |p|^0|. Initially these chains were 
attributed to the chains of tilted vortices, which were 
predicted on the basis of the anisotropic London theory 
fl. Further detailed investigations |l(| indicated that 
this interpretation does not agree with experiment. On 
the other hand the interpretation based on the cross- 
ing lattices qualitatively agrees with the experiment. 
We will calculate the crossing energy of the JV and PV 
stack which allow to obtain a quantitative criterion for 
the transition to the "chain state" . 

First, we consider a structure of an isolated JV. In the 
Meissner state the JV is represented by a region of a 
singular phase distribution []l2 13 1 . The phase difference 



between two central layers changes from to 2n within 
the nonlinear core region. The core size Ajo arL d the 
JV energy e,jo are determined by the balance condition 
between the in-plane phase stiffness energy, ( J/2)(V</) n ) 2 , 
and the Josephson energy, — Ej cos(4> n ,n+i), p3| 
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where 4> n (r) is the phase of the order parameter in the n- 
th layer, (j> n ,n+i = <t>n+i —<f> n — (27rs/$ ) A z is the gauge 
invariant phase difference, A is the London penetration 
depth, 7 is the anisotropy parameter, s is the interlayer 
spacing, J = s$q/7t(47tA 2 ), and Ej — J/{^s) 2 . 

We consider now a JV in the presence of a pancake lat- 
tice with lattice parameter a -C Ajo (see Fig. |l|a) for the 
case of very weak coupling, Aj > A. If we put the JV 
into the PV lattice, its in-plane supercurrents will dis- 
place the PVs from their equilibrium positions. Lattice 
displacements induce extra phase variations, which add 
to the usual regular phase and renormalize the JV struc- 
ture. We will show that at large B z the vortex phase has 
a smaller stiffness as compared to the regular phase and 
gives a dominating contribution to the phase perturba- 
tions. The coarse-grained vortex phase 4> n v(y), created 
by the unidirectional shear lattice displacement in the 
n-th layer, u nx (y), is determined by V y (f> nv — 2im v u nx 
with n v = B z /$>q. The lattice deformations produce the 
phase sweeps between the opposite sides of the deformed 
region 8 y 4> nv = 4> nv (oc) - <j> nv (-aa) related to u nx (y) as 
S y 4>nv = 2nn v J^° oo u nx (y)dy. If we neglect the regular 
phase, these phase sweeps are given by 5 y <j) nv = 7rsign(n) 
in the isolated JV. The vortex phase deformations in the 
JV are large and can not be treated perturbatively. 

To consider a structure of the JV quantitatively we 
calculate the phase stiffness due to the lattice deforma- 
tions. This phase stiffness is mainly determined by the 
shear part of the elastic energy, which for dominating 
magnetic coupling can be written as |I^ , |l5|| 
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with r cut = min(a, A) and r w ui x (0) > s. Using the 
relation between the phase perturbation and the lattice 
displacements, <p v (k) = 2Trin v u/k±, we rewrite the elas- 
tic energy in terms of phase as 
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with the effective phase stiffness J(B z ,k), J(B z ,k) = 
In the range k± < 2/ A and k z > l/r w the 
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i.e., magnetic field reduces the phase stiffness approxi- 
mately by the factor B\/B z . This justifies neglect of the 
regular phase in the region B z ^> B\. The phase stiff- 
ness energy (||) has to be supplemented by the Josephson 
energy. Thermal motion of the PVs induces the fluctuat- 
ing phase 4> n ,n+i and suppresses the effective Josephson 

energy, Ej — > CEj where C = (cos (f> n , n +i*J . Therefore 

the total energy in terms of (f> v is given by 
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It has essentially the same form as in the Meissner state 
except for values of the phase stiffness and Josephson 
energy. Using known structure of an ordinary JV jl2 13 
we obtain for the core size and energy of the JV: 
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The JV energy is reduced by the c-axis field, due to renor- 
malizations of the phase stiffness and Josephson energy. 
The maximum lattice displacement r w and lattice de- 
formation u xy in the vortex core can be estimated as 
r w w aA/Ajo and u xy w aA/A 2 . 

Interactions between the JVs are given by the aniso- 
tropic London theory with the anisotropy ratio 7 e g = 
7y 1 B\j (B Z C). At finite density the JVs form a triangu- 
lar lattice, which is stretched along the layers with the 
ratio 7 G ff, i.e., the JVs form stacks along the c-axis with 
the period c z = y/ /3$o/ '("fe&Bx) and these stacks are sep- 
arated by the distance c y = \/j c s^o/((3B x ). Within the 
anisotropic London theory two states with (3 = 2-\/3 and 
[3 = 2/V3 have the same energy JT(| . These states corre- 
spond to stretching of the regular triangular lattice with 
the closed-packed direction oriented along the layer and 
along the c-axis. The lattice energy Fjl is given by 
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We now compare this energy with the energy change cor- 
responding to the homogeneous tilt of the lattice SF t u t , 
which in magnetically coupled superconductors is deter- 
mined by the tilt modulus cf) at k z = Hpl 
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Comparing Eqs. (|J) and (|J) we see that the crossing lat- 
tices is a favorable state when the tilt angle of the field 
9 B x / B z exceeds the critical value 
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At 9 — 9q the system experiences a first order phase 
transition from the homogeneously tilted lattice to the 
crossing lattices of PVs and JVs. Taking 7 = 500 and 
A = 200nm/ v /l - (T/T c ) 2 at B z = 130 G and T = 70 K 
we obtain #0 w 3°. 

To derive the shift of the melting point by the JV lat- 
tice we write the thermodynamic condition for melting 



F cr (B zm ) + 5F cr (B zm , B x ) — Fi(B z 



5Fi(B zm , B x ), 



where F cr (B z ) and Fi(B z ) are the free energies of the 
crystal and liquid states for field along the c— axis, and 
5F cr (B z ,B x ) and 5Fi(B z ,B x ) are corrections due to 
B x . For small B x we can write B zm (B x ) = B m o + 
SB m (B x ) and expand F cr (B zm ) — Fi(B zm ) with respect 
to 6B rn {B x ), F cr (B zm ) - F[(B zm ) « SB m (B x )AM, 
where AM = Mi — M cr > is the magnetization jump 
at the melting point. This gives 



5B m (B x ) » -AF(B m0 ,B x )/AM 
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where AF(B m0 ,B x ) = 5F cr {B m0 .B x ) - SFi(B m0l B x ). 
Similarly, we can relate the shift of the melting temper- 
ature 5T m {B x ) with the melting entropy jump AS per 
single PV, 5T m {B x ) » -s3> AF(B m0 ,B x )/(B m0 AS). 

SF cr (B m0 , B x ) is determined by the energy of the 
Josephson lattice (@) (we estimated that entropic correc- 
tion to Fjl is relatively weak). In magnetically coupled 
superconductors the melting transition is accompanied 
by strong misalignment of PVs in the different layers || . 
This means that the Josephson coupling in the liquid 
state is strongly suppressed and the influence of the in- 
plane magnetic field in the liquid is much weaker than in 
the crystal, AF(B m0 ,B x ) » F JL (B m0 , B x ) - B^/(8tt) cx 
B x . From Eqs. @ and (0) we obtain that the shift in 
B mz is approximately linear with B x and 
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The linear dependence holds until the JV cores start to 
overlap at B x w 7^/ B\B Z jC. In clean superconduc- 
tors the linear dependence also breaks down at large 
enough B z when the decoupling transition occurs within 
the crystalline phase [[lTj. The melting field of the tilted 
vortex lattice in moderately anisotropic superconductors 
depends quadratically on the in-plane field |l8| in agree- 
ment with experiment JL9| . Therefore the linear depen- 
dence is an indication of the crossing-lattices state. The 
angular dependence of the melting field was measured 
in Ref. H. Fig. |^ shows data from Fig. 4 of this pa- 
per for T= 70K and 80K replotted in the coordinates 



H m cos 9 vs H„ 



H m sin 9. In both cases the 



linear dependence H mz (H m u) is clearly observed. The 
same linear behavior of the melting field has been re- 
ported recently by Ooi et al. M. For T = 70K we have 
B m = 133G, AB = 0.35G, and 8B m /dB x = -0.014. Us- 
ing the same parameters as before we obtain from Eq. (|^) 
an estimate 8B m /dB x w —0.03, in reasonable agreement 
with experiment. 

We now consider the case of small c-axis fields, B z <C 
<f>o/(7 s ) 2 - In this case the JV core contains only one 
row of PVs (see Fig. |l]b). Displacements of PVs do not 
strongly influence the JV structure but create a small 
pinning energy for the PV stack by the JV. Due to this 
energy one can expect at some field a phase transition 
from the weakly deformed triangular lattice to the "chain 
state" represented by regions of good lattice intercepted 
by the vortex chains located on the JV stacks. In the 
"chain state" the degeneracy in the JV lattice [|l6| is 
broken and the state with smaller c z corresponding to 
[3 = 2/V3 has smaller energy. The distance between the 
chains is simply given by the distance between the JV 

stacks c y = yJy/Z^ /{2B x ) @. For B, x w 30G and 
7 = 200 this gives c y ~ 10/im, which roughly agrees with 
the distance between the chains observed in the decora- 
tion experiments P,|lO[|. We now obtain the criterion for 
the transition to the "chain state" . 

First, consider one PV stack crossing one JV located 
between the layers and 1. The crossing energy E x 
appears due to displacements of PVs u n under the action 
of the in-plane currents j„ induced by the JV 
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with j n » gftjfc ( K -i J /'2) 7S > where C n — > 1 at n -> 00. 
For estimates we neglected the weak logarithmic /in- 
dependence of U44(k z ). Minimization of (||) with respect 
to u n gives 
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Numerical calculations give A = J2n°=i ( n C {/ 2 ) ~ 2.1 
and Ai w 3.5. At finite B x each PV stack intersects l/c z 
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JVs per unit length and the interaction energy of the PV 
stack with the JV stack is given by 
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The transition to the "chain state" is expected when the 
energy gain (|To|) to put an extra PV stack on the JV 
stack exceeds the energy loss due to formation of the edge 
interstitial in the Abrikosov lattice Eei at a ^> A [f20| , 
Eei ~ 1.14eo V exp(— a/A). Comparison of the two 
energies gives the following equation for the PV lattice 
spacing at at which the transition takes place 
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The transition field B t — 2$o/v / 3at grows slowly with 
B x . For 7 = 200, A = 220nm, and B x = 30G we obtain 
B t ~ 36G. The transition can be driven by the field or 
by the temperature (due to the T-dependence of A). 

In conclusion, we have shown that in very anisotro- 
pic layered superconductors the vortex ground state in 
tilted fields is given by the crossing lattices of the JVs 
and PV stacks. This ground state explains (i)the linear 
dependence of the c-component of the melting field on 
the in-plane field and (ii)the mixed chains-lattice state 
observed by Bitter decorations at small fields. 
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FIG. 1. JV coexistent with the lattice of PVs. Shaded area 
shows the JV core. Pancake positions in the two central layers 
are shown, a) At large fields the JV core contains a large num- 
ber of pancake rows. The phase field of the JV is determined 
by displacements of the PVs. b) At small fields the JV core 
captures only one pancake row. Due to the interaction en- 
ergy between the PV stacks and JVs, pancake density in this 
row is higher than in neighboring rows (mixed chains-lattice 
state) . 
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FIG. 2. Dependence of the c-component of the melting field 
on the in-plane field for BSCCO at T = 70K and 80K. The 
data are taken from Ref. @|. For both cases a linear de- 
pendence is clearly observed indicating the crossing-lattices 
ground state. 
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